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In double quantum wells electrons experience a Lorentz force oriented perpendicular to the structure 
plane when an electric current is driven perpendicular to the direction of an in-plane magnetic field. 
Consequently, the excess charge is accumulated in one of the wells. The polarization of a bilayer 
electron system and the corresponding Hall voltage are shown to contribute substantially to the 
in-plane conductivity. 

75.50.Pp, 75.30.Ds 



I. INTRODUCTION 

A finite size, in the growth direction, of bilayer electron 
systems confined in double quantum wells leads to strong 
orbital effects of an applied in-plane magnetic field. As a 
result, an oscillation appears in the field dependent den- 
sity of states, corresponding to the depopulation of the 
antibonding subband at a critical field -By = Bc^i, and to 
the splitting of the Fermi see into two separated electron 
sheets at a second critical field i?c,2. The variation of the 
density of states with B\\ is reflected in magnetoresistance 
traces recorded as functions of B\\ . Critical fields B^-u&nd 
Be, 2 can, therefore, be determined experimentalljiiTQ. 

The connection between the magnetoresistance p{B\\) 
and the density of states g{B\\) employed in this papers 
is based on an assumption that the scattering rate is pro- 
portional to the density of states. Unhke g{B\\), which 
is determined exclusively by the electron energy spectra, 
the p{B\\ ) depends also on the electron scattering mecha- 
nism and on the mutual orientation of the electric current 
and the magnetic field. Therefore, the proportionality 
between p{B\\) and g{B\\) is only approximate. 

A conventional approach describing transport in dou- 
ble quantum wells subject to in-plane magnetic fields 
relies on one-electrcwji approximation and semiclassical 
Boltzmann equationcl. This approach is justified in the 
case of "parallel" conductivity a\\ related to magnetore- 
sistance by p\\ — 1/<7||, i.e., when an electric current is 

driven along the magnetic field {j \ \ B^ and no Lorentz 

force acts on electrons. 

Here we present a theoretical scheme which includes 
the electron-electron interactions and allows to evalu- 
ate the "transverse" conductivity a± {p± = l/cr±) cor- 
responding to an in-plane electric current flowing per- 
pendicular to the field direction, j _L B. In this case, 
electrons are driven by a Lorentz force perpendicular to 
the 2D-planes and the electric charge is transferred from 
one well to the other in the current carrying state. The 
Hall contribution to p±, obtained only in our approach, 
is comparable in magnitude to the original p± , and hence 



is important for the analysis of the transport data. 

II. ELECTRONIC STRUCTURE AND 
TRANSPORT COEFFICIENTS 

We consider two coupled, strictly two-dimensional 
(2D) electron layers confined in very narrow potential 
wells at the distance d (see Fig. |l|) . The vector potential, 
A = {zB, 0,0), is used to describe the influence of the 
in-plane magnetic field, B ~ (0,5,0), on the electronic 
structure. 
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FIG. 1. Schematic picture of a double-layer system. Direc- 
tions of applied electric and magnetic fields are indicated. 

Only the lowest bound states of the left and right wells 
are used to construct the model Hamiltonian Hq of the 
systemEl. The Hamiltonian is diagonal in the wavevector 
k = (kx, ky) and takes a matrix form 
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where ^L(R){k) are single- well eigenenergies, 
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and fco = i\e\Bd/2h. The hopping integral t is given 
by a matrix element of the well potential between lowest 
eigenstates of the left and right wells. The diagonaliza- 
tion of the matrix (|l|) yields the bonding and antibonding 
eigenstates. 

We employ a minimal model that accounts for the 
inter-well polarization effect in the Hartree approxima- 
tion. Transport coefficients are obtained using the ran- 
dom phase approximation and Kubo formalism of lin- 
ear response theory. In our geometrical arrangement 
the "parallel" conductivity CT|| relates the current com- 
ponent jy to the external electric field component Ey in 
a standard way where jy = cf\\Ey. The structure of a±^ 
is more complicated as it involves the electron-electron 
interaction contribution. The basic scheme of ourj-calcu- 
lation is outlined below; for more details see Ref. 0. The 
"perpendicular" conductivity is composed of a'"^ , which 
is the response to the applied field component E^^ and 

(2) 

(7]_ corresponding to the "Hall" electric field E^j i.e., 

jx = a'^^ Ex + a^^^ Ez- The "Hall" field appears as a con- 
sequence of the polarization of the sample and is related 
to the inter- well "Hall" potential Uh as E^ — —Un/d- 
It is calculated in a self-consistent way from the non- 
equilibrium excess charge 5Q as a solution to the Poisson 
equation. In our model, 5Q and the potential Uh arc re- 
lated by a particularly simple formula corresponding to 
the parallel plate capacitor. 
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where e denotes the dielectric constant of the barrier. 
The excess charge 5Q is obtained as a sum of components 
(5g(i) and (5Q(2) calculated fr om the Kubo formula as a 
response to the fields E^ and E^, respectively. Writing 



Ex, <5g(2) = £(2)K 
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where e^^^ and e*-^-* are the generalized dielectric func- 
tions, Ez is the solution to the self-consistent equation 
(0) and takes the form 
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This formula relates the self-consistent Hartree field of 
electrons to the x component of the applied external 
electric field E. Using Eq. (||), the perpendicular con- 
ductivity can be written in a form 
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In this equation the first term a^f' corresponds to the 
conventional solution of the Boltzmann equation while 
the second term represents the novel "Hall" correction 
to the conductivity. 



III. SHORT-RANGE SCATTERERS IN THE 
BORN APPROXIMATION 

We have performed numerical calculations of transport 
coefficients for the realistic parameters of a bilayer 2D 
system. Since we have only a limited knowledge of the 
nature of scatterers and their distribution in the sample, 
we assume that electrons are scattered by impurities dis- 
tributed randomly in both left and right wells. To be 
more specific, the scattering on an individual impurity is 
considered to be intra-well (diagonal in the layer index) 
and isotropic in the k space. The concentration of impu- 
rities is assumed very low and the weak scattering on an 
individual impurity is treated in the non-self-consistent 
version of the Born approximation. The strength of the 
scattering on an individual impurity and the impurity 
concentrations are taken as adjustable parameters. 

In the course of calculation we introduced the resol- 
vents in expressions for the dielectric functions and the 
conductivity components. The finite quasiparticle life- 
time and the transport relaxation time are related to the 
imaginary part of the self-energy and to the vertex cor- 
rections to velocity components, respectively. The quan- 
tities result from replacing the resolvents and their prod- 
ucts, characteristic for a given configuration of scatterers, 
by expressions averaged over all possible configurations. 

To proceed further, we adopt additional simplifications 
justified by the very low scattering rate in high-mobility 
samples. Only terms of the lowest order in the scattering 
rate are kept in all expressions. While (T\\, ct^^^ and e'^' 
would diverge in samples without impurities, the leading 
terms, cr^^ and £^'^\ of the expansions obtained as a re- 
sponse to Uh, are finite in this limit. Note that, in spite 

(2) 

of finite (Tj^ , both terms of cr^ are inversely proportional 
to the scattering rate, as the second term is proportional 
to e^^) according to Eq. (||). 

The vertex corrections vanish for the velocity compo- 
nent Vy as expected for randomly distributed short-range 
scatterers. This is not true for Vx- Due to the anisotropy 
of the scattering induced by the in-plane magnetic field, 
the vertex corrections have to be taken into account. 



IV. RESULTS AND DISCUSSION 

The in-plane magnetic field changes qualitatively the 
topology of Fermi contours. At zero field, the Fermi con- 
tours are two concentric circles. The in-plane field shifts 
the centers of Fermi circles in opposite directions in k- 
space and gradually changes the contour shapes. At the 
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critical field -Bc.i, the Fermi line of an antibonding sub- 
band vanishes and the transition from a two-component 
to a one-component system occurs. The Fermi contour 
of a bonding subband splits at the higher critical field 
Sc,2 into two parts and the left and right electron layers 
are decoupled. 
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FIG. 2. Calculated field dependence of the density of states 
g and the magnetoresistance curves py and p±. In p^^* the 
"Hall" contribution is neglected, (a) The concentration of 
impurities is the same in both wells, (b) The ratio of concen- 
trations in two wells is 1:10, i.e. the scattering occurs mostly 
in one of the wells. The schematic shape of the Fermi con- 
tours below, between and above the critical fields Bc,i and 
Be, 2 is also shown 

The calculated magnetoresistance components p|| and 
p± , are presented in Fig. ||, together with the field depen- 
dence of the density of states g. The van Hove singulari- 
ties at critical fields i?c,i and Be, 2 seen on g appear also 
on p|| and p±. Two different distribution of the scattering 
centers are considered, a symmetric case (a) where the 
concentration of impurities is identical in both wells, and 
an asymmetric case (b) where the prevailing amount of 
impurities is in one of the wells and the ratio of their con- 
centrations is 1:10. Three regions of the magnetic fields 
should be discussed separately: below -Bc,i, between Bc.i 



and i?c,2 and above i?c,2- 

In Fig. Ija, the curves are strikingly similar below i?c,i, 

except of p±. The components p|| and p^^\ which can 
be obtained by solving the Boltzmann equation, exactly 
coincide with g for our simplified model. It is the Hall 
contribution which qualitatively changes the field depen- 
dence of p± to negative magnetoresistance in this re- 
gion. Between Bd and i?c,2 only the bonding subband 
is occupied and the difference between p|| and p± origi- 
nates from the Fermi contour anisotropy and the related 
anisotropy of the effective mass and the relaxation time. 
The anisotropy of the transport relaxation time (the ver- 
tex correction to v^) plays a role also above i?c,2- 

Fig. ^ shows the magnetoresistance curves calculated 
for the asymmetric distribution of scattering centers. 
The dramatic reduction of all magnetoresistance com- 
ponents with increasing magnetic field is related to the 
in-plane field suppression of the coupling between wells. 
The tunneling between wells is reduced with increasing 
magnetic field and the conductivity of the well with lower 
concentration of impurities prevails. In this case the 
"Hall" correction does not change the field dependence 
qualitatively but enhances the difference between p|| and 
p± substantially. 

To conclude, we have shown that the charge redistribu- 
tion within a double-layer system gives rise to the Hall- 
like contribution to the perpendicular magnetoresistance 
component p± which is so strong that it cannot be omit- 
ted in the analysis of bilayer 2D system magnetotrans- 
port. 
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